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Within the context of finding the initial conditions of the universe we consider gravitational in- 
stantons falling into the Bianchi IX classification. That is, a Euclidean four-manifold with a metric 
that satisfies Einstein's equations with an induced metric on S 3 submanifolds that is homogeneous 
but anisotropic. As well as finding regular solutions to the field equations with a tunnelling scalar 
field, we also look at the case of singular instantons with a view to applying the results to generic 
potentials. The study is in agreement with the prejudice that instantons with higher symmetry have 
a lower Euclidean action, even when we consider the singular class of solutions. It is also found that 
the Euclidean action can diverge for simple potentials, showing that the Hawking Turok instanton 
had finite action owing to its symmetry. 
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I. INTRODUCTION 

Questions of initial conditions in the universe have 
been around for many years, but only now with the 
promise of detailed data from the cosmic microwave back- 
ground can we expect significant experimental input. 
The proposal which this work is concerned with is that 
of the universe 'tunnelling from nothing'. That such a 
situation could be realised was noticed by Vilenkin El, 
where he pointed out that the compact nature of the 
de Sitter instanton could be interpreted as the universe 
coming from nothing into a de Sitter phase via the de 
Sitter instanton. In a similar vein Hartle and Hawking 
P| attempted to describe the initial conditions of our uni- 
verse by proposing that the Lorentzian metric of space- 
time be rounded off with a Euclidean manifold. The 
two approaches are similar in philosophy but different in 
details, leading to different probability measures for the 
initial conditions of the universe j| . 

Studying instantons is at some level independent of 
these considerations. What one is looking for is a Eu- 
clidean solution of the field equations which allows a 
Lorentzian solution to be glued smoothly onto it, the 
question of probability measures comes later. As such, 
we start by just trying to find solutions to general rela- 
tivity on a manifold with Riemannian metric, a daunting 
task as there are an infinite number of degrees of free- 
dom in the system. To aid us we impose some symmetry 
requirements that leave us with a tractable problem, this 
is not entirely ad hoc as the analogous system of a field 
tunnelling in flat space shows. Coleman et al H showed 
that imposing SO (4) symmetry on a bubble at zero tem- 
perature was sensible in that it had the lowest action and 
was therefore more likely to nucleate. In fact it was later 
shown that it was more than just sensible, because all 
other solutions to such a system are singular H . 

The idea of singular solutions has been brought to the 
fore by Hawking and Turok jq]; there it was pointed out 



that whereas singular solutions in flat space means in- 
finite action, this is not the case in gravity. Including 
gravity means that when the solution starts to blow up 
the manifold at least has a chance to close up and make 
the volume of the instanton finite, allowing a non diver- 
gent action. Whether such solutions are to be allowed 
into the path integral has been hotly debated, here we 
take the opinion that we should really let the action de- 
cide what is allowed rather than introducing some notion 
of singularity by hand to disallow these instantons. We 
shall find in fact that this is more restricting than may 
be expected. It was asserted in g that generic potentials 
lead to finite action for solutions with SO (4) symmetry 
(although this was later shown to be true only for poly- 
nomial and certain exponential potentials M). We will 
find here that even with polynomial potentials the sym- 
metry of Bianchi IX is not always enough to keep the 
action regular, unlike imposing SO (4) symmetry. 



II. EINSTEIN EQUATIONS. 

In this paper we consider the field theory of general rel- 
ativity coupled to a scalar field according to the following 
action, 
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boundary term, (1) 



where the boundary terms are designed to cancel the sec- 
ond derivatives of the metric, allowing a sensible meaning 
to the variational principle on a boundary J8|], ||. The 
gravitational coupling is K — 8nG, which we rescale to 
unity, and r\ is the volume four-form of the manifold. The 
metric is taken to be positive definite and of the form, 

ds 2 = dr 2 + a(r)V) 2 + b(r)V) 2 + c(r) 2 (a 3 ) 2 . (2) 

The o % are the left invariant one forms of SU(2) and so 
satisfy da 1 = — ie 1J <r J A <r k , making the solution fall into 



the Bianchi IX class. Using Cartan's structure equations 
for a torsion free connection it is a nice exercise to find 
the curvature two forms, leading to the following field 
equations, 
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Here (3{t) = a(T)b{T)c{r) is the volume of the surface at 
constant r and ' denotes differentiation with respect to 
t. To get the equations of motion for b(r) and c(t) we 
use the Z3 symmetry of the equations and just cyclically 
permute a, b and c in (0). 



III. THE REGULAR SOLUTIONS. 

We know that for a potential with a false vacuum there 
is the Coleman-De Luccia solution pOf where we take the 
functions a(r), b(r) and c(t) to be equal, the isotropic 
limit of Bianchi IX. When we want a regular solution 
we require it be regular everywhere, in particular at the 
origin. In the coordinates we are using one finds that 
the Coleman-De Luccia solution has the scale factor be- 
having as a(r — > 0) — > r/2 giving the instanton a metric 
which locally is like an isotropic metric on R . This how- 
ever is not the only way to smoothly end an instanton. 
The Coleman-De Luccia solution ends on what Gibbons 
and Hawking |LlJ called a nut, where all the scale factors 
vanish linearly as \t. There is also the case of closing 
the instanton off with a bolt, in this case where one of 
the scale factors vanishes linearly as |r the other two 
become equal and constant. One could also consider the 
situation where the vanishing scale factor does so with 
unit proportionality to r, this changes the group orbits 
of the homogeneous submanifolds to SO (3) rather than 
SU(2). Now we consider the regular bolt solution with 
SU(2) orbits on the homogeneous submanifolds, calling 
the point at r = the south pole. At the bolt then we 



impose that c varies as ^r and a, b are equal and con- 
stant, with the equations of motion showing that they 
remain equal. 

A nice geometrical picture exists for these instantons 
in terms of expanding and contracting three spheres. 
The constant r slices are seen to be homogeneous three 
spheres, whose anisotropy is dictated by the scale fac- 
tor at that time. The nut-nut solution is viewed as an 
isotropic three sphere expanding from zero size at the 
south pole and contracting back to zero at the north 
pole. Looking at the bolt-bolt case we see that the ra- 
tio of the unequal scale factors represents the anisotropy 
of the S 3 slices, which diverges at the poles. This limit 
of anisotropy is where the three sphere, viewed as an S 1 
bundle over S 2 , has fibres of zero length, giving the two 
sphere limit of a squashed three sphere. The picture we 
have then is that at the south pole we have a maximally 
squashed S 3 i.e. S 2 , becoming less squashed as we move 
off the pole, finally becoming a two sphere at the north 
pole. 

With the three volume [3 increasing linearly from zero 
at the south pole bolt we may infer from the Euclidean 
Raychaudhuri equation [[l2| that for some finite r > [3 
will again go to zero, assuming that the scalar potential 
is positive semi-definite. The point at which the three 
volume goes back to zero we shall call the north pole. The 
manner in which it closes, i.e. singular or not, depends 
on the value of the constant that a(0) (= 6(0)) is chosen 
to take, as well as on the value </>(0). We can see that 0(0) 
is not fixed by the constraint equation (]4]) by noting that 
near the bolt a — ► x + C 7 " 2 1 with \ an d £ arbitrary. The 
constraint equation then relates 0(0) to these arbitrary 
constants. Given that the 3 volume, /3, vanishes and we 
are looking for regular 0(r) then (gj) shows us that 0' = 
at the north pole, as (3' / (3 diverges there. 

A similar argument to the undershot-overshoot argu- 
ment of |l|] suggests that values of 0(0) and a(0) exist 
such that the instanton ends on a bolt at the north pole, 
if we have a potential with a false vacuum at = and 
true vacuum for some positive 0. For a given a(0) we 
may take 0(0) to be close to the vacuum value, in which 
case will diverge to negative values as we approach the 
north pole. If we start further from the vacuum value 
then will have turned around by the time the damp- 
ing force f3' / f3 of (ra) becomes anti damping, driving 
to positive infinity. There is thus a value between these 
which make the scalar field constant at the north pole. 
We then find however that a(0) has not been fixed yet. 
At the north pole the functions are now regular, however 
it will not in general be a bolt, with c(t — » tn) going like 
7(t — tn) and 7 being different from a half. We now use 
the freedom in the value of the constant that a(0) takes 
to find the solution with J = h- 

To make these arguments concrete we took the poten- 
tial 

V(0) = (0 2 -l) 2 (0 2 +a), (6) 

which exhibits a false vacuum at the origin for < a < ^ > 



we took a = 0.1. A regular solution using this potential 
has been found numerically and is displayed below in 
Fig. 1. We see from the profile of <\>{r) that the scalar 
field interpolates between the false and true vacua, as in 
the Coleman-De Luccia case. However in that case one 
is able to find a section of the instanton which has zero 
extrinsic curvature, allowing one to use this surface to 
attach a Lorentzian metric to it corresponding to an open 
universe. Here however we do not expect to find such a 
three-surface owing to the anisotropy of the homogeneous 
submanifolds. We can however find a solution which has 
an inversion isometry, implying that the stable surface of 
this symmetry has zero extrinsic curvature J[4f. Such a 
solution is shown in Fig. 2 where we see that the profiles 
are symmetric about the mid point between the north 
and south poles. Such an instanton nucleates a closed 
universe with anisotropic spatial sections. 
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FIG. 1. profiles of the functions involved in a regular bolt 
instanton. 




IV. THE 'CONICAL' SOLUTIONS 

The suggestion that singular solutions are valid allows 
us to consider the case where the metric functions are reg- 
ular but the curvature diverges, reminiscent of the cur- 
vature blowing up at the end of a cone. This option was 
not available in the isotropic case, as the constraint (Q) 
fixed the gradient of the metric functions. The freedom 
of the bolt type boundary conditions allows us to satisfy 
the constraint with c(t) varying linearly with any gradi- 
ent and a(r) (= b(r)) being an arbitrary constant. Here 
we consider such a case, imposing the condition that the 
instanton have a parity isometry such that we may at- 
tach a real Lorentzian metric to the stationary surface. 
On this surface the scale factors a and c take some con- 
stant value and their ratio, ( = a(T sym ) j c(T sym ) , gives a 
measure of the anisotropy, such that the isotropic case 
has £ equal to unity. We then expect that varying £ 
will allow us to move continuously from the bolt-bolt in- 
stanton to the Coleman-De Luccia solution with parity 
isometry. This is indeed the case and is illustrated be- 
low in Figs. 3,4. The upper set of curves in Fig. 3 are 
the a(r) scale factors as we vary the anisotropy measure, 
C, whilst the lower set are the c(t). We see that one 
is able to interpolate between the two regular solutions 
(i.e. the Coleman-De Luccia case with a(r) = c(t) and 
the bolt-bolt of Fig. 2) by using the solutions which have 
a conical singularity. To make the comparison clearer in 
this figure we have scaled all the solutions to have the 
same size, whereas in actual fact the instantons become 
larger as we approach the isotropic limit. This informa- 
tion is given in Fig. 4, where we can see how the size, 
tn, and the Euclidean action (including boundary terms) 
varies with anisotropy. One can see that the action is 
lowest for the isotropic case, confirming naive expecta- 
tions. We also see that the bolt-bolt instanton has no 
privileged position, the anisotropy of the bolt-bolt case 
is around £ = 2.4 at which point the action shows no 
special behaviour. We also note here that the dotted 
line corresponds to the, separately calculated, nun-nut 
(Coleman-De Luccia) case. 
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FIG. 2. profiles of the functions involved in a regular bolt 
instanton with a parity isometry. 
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FIG. 3. Interpolating between a nut-nut and bolt-bolt in- 
stanton. The upper set being the a scale factor and the lower 
being the c. 







-2 



"N 



action 



1 1.5 2 2.5 3 

anisotroov 

FIG. 4. Action and size of the instantons which interpolate 
between the nut-nut and bolt-bolt solution. 



V. THE SINGULAR SOLUTIONS. 

We now move on to consider the solutions which pos- 
sess a singularity on the manifold different from the con- 
ical type discussed above. The type of singularity we are 
considering here is the analogue of the singularity stud- 
ied in the Hawking Turok paper H. This singularity is 
different in that the curvature here diverges as l/(r— tn) 
whilst the conical singularity is more like a delta function 
at the poles. 

As we are not looking for regular solutions we drop the 



tunnelling form of the potential used before and work 
with a simple quadratic potential. First we look for 
solutions which contain a reflection isometry. This is 
achieved by setting initial conditions at some point such 
that the fields a, b, c and (f> all have vanishing gradient 
at that point, this allows a closed anisotropic universe to 
be grafted on to this surface. 

We will be interested in finding the Euclidean action 
for such instantons and as such we need to take care 
of the boundary terms mentioned in (fil). On our solu- 
tion the singularity must be excised from the manifold, 
this means taking a small region around the singular- 
ity and discarding it, naturally introducing a boundary. 
The boundary term is then the contribution this bound- 
ary makes as the region is made arbitrarily small. We 
find that the Riemann curvature for the metric (0) is, 
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giving the boundary term of ([jj) as, 



boundary term = - (/3'(r -* t s ) - 0'(t -> t n )) V 3 (8) 
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In order to calculate with such ill behaved solutions one 
needs to be confident of the numerics. To be sure that 
the correct answers are at least well modelled we ran a 
number of simulations with different integrators. The 
results presented below are those calculated using a Bul- 
lirsch Stoer method E3|, which was tested against a 
Runge Kutta routine also from jlj| and the integrators 
found in mathematica and maple. All gave the same re- 
sults to fractions of a percent. What one finds in do- 
ing this is that this simple potential does not always 
give finite action, as it did for the 0(4) invariant solu- 
tion. Following []5| we assume that near the singular- 
ity the form of the potential becomes irrelevant, a fact 
backed up by numerical solutions with polynomial poten- 
tials. The equation for <f> then shows that near the north 
pole 4>'(t — > tn) — > 1//3. Taking the constraint equation 
(Q) assuming a form of a — > (t^ — t) a , b — » (tjv — t) b , 
c — > {tn — T ) C \ where A, B and C are positive, we then 
drop the terms which are products of scale factors. Us- 
ing the fact that <f>' is varying inversely as the three vol- 
ume one finds A + B + C = 1 , showing that dropping the 
products of scale factors was consistent. The three vol- 
ume, /3, therefore goes to zero linearly, implying that the 
scalar field diverges logarithmically, which is slow enough 
that the action can converge. 



When looking at the action as a function of initial val- 
ues for a, b and c we found the very striking result that 
unless two of the scale factors were equal the boundary 
term would diverge. Inspection of the numerical solu- 
tion revealed that for such solutions not all of the scale 
factors went to zero at the singularity, this effect is com- 
ing about because it is the scale factor products, rather 
than <f)' which dominates in the singularity. In fact such 
behaviour can also exist when two of the scale factors 
are equal, say a and b. It is found that the boundary 
term diverges unless c(t = 0) < a{r = 0). It is interest- 
ing to note that this biaxial behaviour where two scale 
factors are equal also appeared for the regular solutions, 
but for different reasons. The regular instantons were re- 
quired to be biaxial, extending the isometry from SU(2) 
to U(2), so that the instanton could be closed off at a 
regular bolt, see also |1J|. The picture we then arrive 
at is given below in Fig. 5, which is a plot of the total 
action, when it is finite, for a set of singular instantons. 
The figure shows a set of contours indicating that the 
minimum action occurs for the isotropic instanton with 
o(t = 0) =6(t = 0) = c(t = 0). 
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c(x=0) 
FIG. 5. Total action for the singular instanton with reflec- 
tion isometry and b(r) = a(r), only the finite action solutions 
are shown. It is seen that the minimum action occurs along 
the line of isotropy. 

Although this figure only covers the b(r) = a(r) plane, 
it appears from the numerical results that it contains all 
the relevant information. The Z3 permutation symmetry 
of the scale factors mean that this plane is just the same 
as the 6(t) = c(t) and a(r) = c(t) planes, with the action 
diverging off these planes. We then note that the minima 
of the action is to be found along the isotropy axis, that 
is to say the metric which has the 0(4) isometry and is 
the gondola instanton studied in m& . 



VI. CONCLUSIONS 

We have discussed a generalisation of the SO (4) sym- 
metric instantons, keeping the homogeneity assumption 
but relinquishing isotropy. In particular, our solutions 
were of the Bianchi IX type. It was found that the sin- 
gular instantons in this class which allowed a reflection 
isometry, leading to the initial conditions of a closed uni- 
verse, had a minima of the action when the instanton was 
isotropic. Moreover, the probability distribution had a 
Z3 symmetry around the axis of isotropy so the average 
instanton would lie on this axis (both mean and mode). 

An in depth search of types of potential has not been 
carried out here, concentrating only on the simplest of 
potentials. Indeed, the results of this study did not reveal 
any reason to expect other potentials to deviate from the 
behaviour found here. This applies to both the regular 
solutions where the potential has a false vacuum and the 
singular solutions for generic potentials. One of the main 
results is that even for the singular instantons isotropy 
would seem to be preferred over anisotropy. The study 
also sheds light on the issue of whether or not to allow sin- 
gular instantons, in that the action can in fact be a good 
discriminator of solutions, which may mean that remov- 
ing solutions by hand is not necessary. Namely, we saw 
that it was the high degree of symmetry in the Hawking 
Turok solution that meant the criterion of finite action 
would keep all solutions with polynomial potentials. 
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